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ABSTRACT 


A simple theoretical relation between a one-dimensional sur- 
face wave and the amplitude fluctuations of an acoustic standing 
wave in a cavity with idealized boundaries is presented. An ex- 
periment was designed to investigate the applicability of the 
theory to a more complicated configuration. A 237 cm x 117 cm 
x 9 cm water-filled cavity with pressure-release walls was excited 
in a normal mode by a source placed in one corner. The resulting 
acoustic signal was monitored at various locations throughout the 
cavity with a small probe receiver. Wind-driven surface waves 
were generated by a multi-fan blower placed at one end of the 
cavity. The relationship between the resulting fluctuation of the 
acoustic Signal and the surface fluctuation above the receiver was 
investigated. The results are consistent with the simple theory 
at low surface-wave frequencies. Anomalous behavior was observed 
for the higher frequencies of the surface wave spectrum. Possible 


mechanisms for this inconsistency are postulated and discussed. 
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I. INTRODUCTION 


In A.B. Wood's model study of shallow water propagation Hilly a 
sound amplitude fluctuation, synchronized with the passage of sur- 
face waves, 1S reported. Wood suggested the cause was an appre- 
ciable change in water depth in the vicinity of the source and/or 
specular scatter from the surface. The reflection and scattering 
of sound from uneven surfaces has been given considerable attention 
in the past twenty years [2]. Most treatments in the open litera- 
mee have involved ray theoretic approaches, treating the reflect= 
ion and scattering of plane waves from irregular surfaces, and 
have assumed no interaction with a second boundary. When a second 
boundary is present, the possibility of ducted propagation arises. 
A practical mathematical tool for investigation of surface irregu- 
larities affecting this type of propagation would be normal mode 
theory. Normal mode approaches [3,4] have been rare. 

Several experimental investigations involving normal mode pro- 
pagation have been conducted. J.A. Scrimger [5] reported corre- 
lation between surface-wave heights and sound-amplitude fluctuat- 
ions both in magnitude and in frequency content. 0O.S. Tonakanov [6] 
found a linear relationship between surface-wave heights and per- 
centage fluctuation in sound amplitude and phase. 

The present work describes an experimental method for the in- 
vestigation of the effect of surface EGA Aion on sound-fluctuat- 


10n in guided mode propagation. Sound-amplitude fluctuations, 








caused by a wind-driven surface, were measured and compared to 
the surface irregularity in, the vicinity of the 7eceiver eine 
surface was treated as a Fourier superposition of varicus sinu- 
soidal waves, and the experimental results at each frequency were 
compared to a theoretical prediction for sound-amplitude fluctu- 


ations based on a monofrequency sinusoidal wave surface. 
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ieee LHe Oy 


A normal mode solution to the wave equation for a rectangular 


enclosure with smooth boundaries (Figure 1) can be expressed as 


$ 


I 


cos k,x cos k,y sin k,2 cos Wt (Zale) 


Cimeularetrequency Of excitation 


I 


where WwW 


; cthoee 
DpeOpegat1On Constant tor the 1—— divection: 


i] 


Ke 
aL 


If the boundaries perpendicular to the z axis are pressure- 


release and are located at z = O and z = Lo then k36) = nN, is 
required. In all that follows, it 1S assumed that the system is 
excited at frequencies allowing only ny = 1] so that the eigenvalue 


k., is as small as possible. 
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Figure 1. Theoretical Cavity for Guided Mode 
Propagation 
If the upper boundary is disturbed, the boundaries in the z 
@arection then occur at 


z= 0 


Z Le spe ices i ae (2a) 
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where €F(x,y,t) describes the change in Gepth of the cavity re- 
sulting from the surface fluctuation. 


If SF(x,y,t) as small, then a Solwmtiloneo! sinesier) 
o = + € , 
1+ ef, (2-3) 


is sought where 2, is the solution for smooth boundaries shown 
above and cB, expresses a perturbation due to the altered boundary 
ae Zz = £ 
a 

When € is sufficiently small, the following boundary conditions 
on B,, result L7- 
0G, 
a ee a) (Ze) 


Ze Pas, 
Z Z 


Q 
I 


z=0 


and 0, must satisfy the wave equation 
a 
L] 9, =i, 


If the surface fluctuation is a monofregquency wave propa- 
gating across the surface z = Le in the x direction, F({x,y,t) may 
be written 


F(x,y,t) = ss cos (y¥x-Mt) (225) 


where Y = surface wave propagation constant 


{2 = circular frequency of the surface wave. 
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The Solute for fie. under these conditions 1S 


ae S iter 
o, ai Bee lis Koy ST cosl (k, +7) x~ (wd) t] + (235G)) 
5 2 
Sin 2 
—_— - Ut? 
= ra cosl (k, yay saa )t J+ 
sin oe 
SeL) cosl (k, #7) x+(W-Q) t} + 
eZ 
sin ne 
ae cosl (k, ~Y)x-(@) t] 
ui 
2 THOe 2 2 
where Da = aa) - (Kk, +Y) - Ky (2.7) 
2 ant? e : 2 K 2 
an a ( =e) 7 ( 17”) iz 2 
2 UW ‘ ae 2 . 2 
P39 a ( C ) ( ge a Z 
2D w-0) : Z 2 
an = ( Cc ) ~ (kk, -¥) a Ky 
If Q<< w, then the approximation 
2a 2 Z 2 
Pe oe ag) ak, (25) 
a > oe 2 
ae) = sek = k ~(k, -¥) -k, 


can be made and the complete solution simplifies to 


@ = cos Wt cos kyl cos k x Sin ke2 +s ek Acos(QtH)) J (2eo9) 


3 


where 









: Z a : | : 
sin Ib 2 sin Woz Sin 52 sin Ib 2 
sin tb_% Sin 2% bios Sin w_L Sin 4 as ae 
ZZ Zee TZ 2 2 


AF. 








and 








‘sin M42 sin 2 
Scat + { —p-——S— |] si 3 

ania sin(k, ys eT) sin(k, y)x 

lz 2 
Tam) = - : a 

sin b,2 sin [2 

Santi b. poo tai) )* | eamerneey | COS amas 
lz Ze 


Notice that it is the second term in the brackets of Eq. 2.9 that 
describes the fluctuation in the acoustic amplitude and that this 
term depends only on x and z, For the subsequent analysis Eq. 2.9 


mel be written in the form 


¢ = cos Wt cos Ky [cos kx Sin SS + @ cos (Qt+w) ] (2 22n0;) 
where kf = 
rez 
a/— =(T/2yA. 


When kx =m (an aniinode of the standing wave in xX), 


a/E becomes 











: 2 ; . : 
Sin Wz sin jz sin Wz Sin pz 
a/e = See Ac Sate. tee! S i he i Sin wU£ (2.11) 
1 z 1 2 OD: : Be 
At z = hie, this expression further reduces to 
” J. 1 
Be 7a cos pk * Cos bad ceca) 
Le 22 : 
Z 2 
and 
®€ = cos Wt cos ky Lita cos(Qt + )]. (2283) 
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Predicted Values of a/c 
Vs 
Surface Wave Frequency (2/27) 
{at an Antinode of the Standing Wave) 
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Figure 2 
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Predicted Values of a/e 
VS 
surface Wave Frequency (2/2n) 


(at a Node of the Standing Wave) 
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Figure 3 
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On the other hand, at a node where kx =TN/2 ana for 


x = Lo/2, Eq. 2.12 becomes 





it He 
oy i es (2.14) 
£ Me 
cos pk) cos Lp : 
fad Zz 
Baad 468 )6= a cos WE Cos ko cos(Qt+ a). (22185) 


Plots of a/€ as functions of surface wave frequency for the 
two cases are presented in Figures 2 and 3. These functions dis- 
play strong dependence on BOEn acoustic and surface-wave propa- 
gation constants. 


Notice that when Y = k there 1S A= resonance ~ciiccre. pine 


1? 
perturbation analysis may fail in this region, because evr may no 
longer be small compared to 9): If bulk absorption of the medium 
is included, a/€ no longer becomes infinite, but is still extremely 


Maroe. For y < k the surface fluctuation tends to behave as a 


ile 
rather uniform change in depth of the cavity, and the amplitude 
mittetuation behaves as if the standing wave in the z direction 


were squeezed and then released. When Y 7 k the amplitude 


1? 

fluctuation appears to react to the average fluctuation over many 

Surface wavelengths, and approaches the steady-state solution ,- 
The relation between frequency and propagation constant is 


W/k = c (2.16) 


for the acoustic wave and 


O/y V (22) 


for the surface wave where c and V are the respective phase speeds. 


ile, 








The surface-wave speed V is dispersive and the governing 


equation [8] is 





2 Wee 
vo=(d/y + A ) tanh (Ye) (21a) 
: . | ‘ Z 
where Gg = acceleration resulting from gravity 980 cm/sec 
Oo = density 
. ee 
Te = surface tension 75 dynes/cm 


Because of this, Y varies with 92/27 as shown in Figure 4. 
In this experiment, {2 is of magnitude 10 rad/sec and 


i; 


oe 90 x 10? vad/see. “Therefore, VV ~ ko anda —9o0n. 108. an 


1 
mite teason, the perturbation term in Eq. 2.13 appears as a low- 
frequency modulation of the acoustic signal, and @ becomes a 
measure of the modulation amplitude. The great inequality in 
acoustic and surface wave speeds means that the surface can be 
treated as effectively "frozen" with respect to the acoustic 
processes. 

The solution for eG, in terms of €F(x,y,t) obeys super- 
position. Thus, if the surface were a spectrum of §2, the solution 
ieOr EB, would be superposition of the solutions for each component 
Gene according to Eq. 2.13. 

In particular, for a narrow band of surface wavelengths, a 
localized disturbance on the surface leads to a local disturbance 
in the acoustic standing wave. Scrimger [5] observed that the 
fluctuations in the acoustic signal were strongest when generated 
surface waves were in the vicinity of the source or receiver. This 


localization of effect is assumed to apply in this experiment. 
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IIIT. PROCEDURE AND ANALYSIS 


Peer EXPE RKIMENE 

The goal of the experimental investigation was to achieve 
quided-mode propagation in a shallow water layer, disturb one 
surface with a measurable perturbation, and observe the resultant 
fluctuations in the acoustic pressure amplitude. 

In the original concept of the experiment design, one end of a 
shallow tank with pressure-release boundaries was to be lined with 
sound absorbing materialandafan manifold for generating wind- 
driven "seas'' was to be suspended over the other end. By using a 
Peohly directive source aimed into the anechoic area, a cylindri- 
cally spreading acoustic traveling wave, which excited the n3 = ] 
family of modes, could be transmitted. This approach would have 
yielded the following advantages: 

1. Elimination of effects resulting from reflections 
at the vertical sides of the tank. | 

2. A Simple traveling wave configuration corresponding to 
propagation in the lowest mode of a cylindrically 
diverging wave. 

bar limitatilon of the ensenified surface to a xregion 
excluding that in the vicinity of the fans used to 


generate the surface fluctuations. 
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It was impossible within the available time to construct a 
sufficiently anechoic tank in the frequency range of interest, so 
the method had to be abandoned. The alternative was to place a 
cylindrically symmetrical source in a corner of the tank and ex- 
cite a normal mode of the bounded cavity. This method raised the 
immediate problem of the effects of the highly agitated surface 
directly under the fans..(See Figure 5.) (The surface under the 
fans was beaten about by the air strean generated by the fans 
directed at the water surface.) 

The theory suggests (as discussed earlier) that it may be. 
possible to assume that for small €F(x,y,t), the cause of the 
acoustic amplitude fluctuation observed at the receiver is the 
Eeetace in the immediate vicinity. Therefore, acoustic noise fron 
the turbulence under the fans could be treated as a "distant 
source" combining incoherently with the effects generated by the 
local surface. 

Excitation of the ge 1 family of modes was a primary con- 
Sideration; the complications coming from model interference were 
to be avoided. To ensure efficient excitation of these modes, a 
source active over a height greater than the depth of the cavity 
was chosen. This insured that, as the depth at the source changed 
with passing waves, source strength exciting the desired normal 
mode would be smoothly changing and closely related to the local 
surface. Since sea state at the source was essentially uncor- 
related with that at the receiver in the cavity, and since the 


propagation speed of the acoustic energy far exceeded the 
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speed of the surface waves, the effect was the equivalent of 
simply changing the amplitude of 2, - Equation. 24! 3.nlatethen 


be modified £0 


6 = Ll+a. cos(8lt +  _)]cos Wt cos kyl 1+a cos(Qt+wv)] (3a) 
where a is the amplitude fluctuation at the receive resulting from 
the changing conditions at the source, and De and W are assumed in- 


coherent. For small a. and a this becomes 


6 = cos Wt cos kyl 1+ Na.” + ae cos(2t++ vty]. C302) 


In a Similar manner, the "distant acoustic sources" contributed 
by the strong surface-fluctuations in the neighborhood of the fans 
combine incoherently at the receiver with the acoustic amplitude 
fluctuations resulting from the local surface state. Inclusion 


of these "distant sources" modifies Eq. 3.2 to 


6 = cos Wt cos k,yll+ | 074 a” + a. cos(Qt + wy] (3ion) 


where a, = the emplitude fluctuation at the receiver resulting 116m 
the distant sources. 


Thus, the sound-amplitude fluctuations resulting from the local 


surface state must be obtained from the measured fluctuations and 


\ Z Pa Zz 
a. = a oe as + a. (3.4) 





a = Va - a, - a : (ao) 


Extraction of the desired values of a and then € from the ex- 


perimental data is accomplished by the following steps: 
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1. Measure the statistics of the sea surface to obtain 
€ as a function of surface wave frequency (2); 

2. Measure the acoustic amplitude fluctuations (a.) at 
the various positions in the tank; 

3. Determine the values of Qs and a, as FUuNnecteronsweL 
surface wave frequency; 


ie Compute ¢ by Equa 3.5 and vdiviide at) by Cua tieaen 


frequency of interest. 


Ee DESCRIPTION 

The experiment was conducted ina shallow water tank whose 
dimensions are shown in Figure 6. The tank was constructed of 
3/4 inch marine plywood set on 2 inch x 6 inch supports to pre- 
vent sagging. Water depth could be maintained uniform to within 
fee mm throughout the tank, 

7.3/3 inch no enee envelope, set in the tank and held in 
place by hydrostatic pressure, served as a pressure release 
boundary. The Fate coefficient of the neoprene was measured 
by terminating an acoustic waveguide with a sample of the material 
and conducting Berenice tests. Walves Of around 0205 106, the 
reflection coefficient were obtained. The boundary condition in 
the tank was confirmed by probing the sound field and verifying 


that the pressure extrapolated to zero, within experimental accu- 


racy, as the probe was brought near a boundary. 





The neoprene rubber used is foam containing nitrogen and is 
the type used in suits for SCUBA divers. It is commonly known as 
"wetsuit" material. 


2k 








a en ee eee 


ky en 1Op 


Inside Dimensions 
Fiberglass 
Coating 


3/4 inch Marine 
Plywood 


LLLLLLLL 


2 3/8 inch Neoprene 
~~ ¢—— Rubber 


WO 


Tank Construction 
- Figure 6 
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A "wind-driven sea'' was created by five equally~spaced 
squirrel-cage blowers in a manifold at one end of the tank. The 
manifold was suspended over the tank to prevent coupling of 
mechanical vibration to the tank (Figure 7). Each fan had a rated 
capacity of 100 £t7/min. Various seas could be generated by varying 
ae angle of attack of the fan output, the height of the manzieola, 
orifice size, and the number of fans (Figure 8). Greater sea 
heights were obtained by placing a 16 inch wide cover over the 
water adjacent to the fan Aerekaouta in order to hold the wind down 
on the surface (Figure 5). 

The wave heights were measured by an immersion probe whose 
impedance varied inversely with its submerged length. This wave 
probe consisted of two submerged electrodes with a 10 kHz cur- 
rent flowing between them. The passing surface waves modulated 
the current amplitude by changing the amount of submerged elec- 
trode. This produced a current whose envelope amplitude was a 
measure of the passing wave height. After amplification, the 
information was retrieved by filtering and detecting the amplitude 
envelope of the current. 

Wave probe construction is shown in Figure 9. The positive 
electrode was made from 20-mil nikomel wire soldered to 1/16 inch 
brazing rod. The gound return was 1/16 inch brazing rod coated 
with neoprene paint over its vertical length. The last 4 cm were 
kept bare and bent at a right angle to the vertical. The elec- 
trodes were soldered to standard banana plugs and held in place 


by a connector (Pomona Electronics Connector model #1681). 
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Tank with Fan Manifold in Place 


Figure 7 
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Wind Generator (Fan Manifold) Operation 
Figure 8 
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The Wave Probe 
Figure 9 


END) 





The detection and filtering network (see Figure 10) con- 
sisted of a switching diode for half-wave rectification, two 
sucessive low-pass filters with cutoff frequencies of 159.5 Hz, 
and a blocking capacitor. The 25 UF capacitor with the 50 kf 
input impedance of the following amplifier formed a high-pass 
filter with 0.126 Hz cutoff. The network was tested with an 
unmodulated 10 kHz tone which ideally would generate zero DC and 
AC outputs. The AC signal which leaked through the filter amounted 
momo more than 0.05 percent oe the input amplitude, and the DC 
appearing at the output was equally ImSiondf teat. 

The relation between submerged length of the wave probe and — 
current flow was non-linear, the effect being more or less severe 
depending on the positioning of the probe with respect to the 
cavity boundaries. (If the electrodes were placed too near the 
bottom, the non-linearity was aggravated.) By bending the ground 
return and adjusting the probe current level, an adequately linear 
operating range was established for the wave heights expected. 

Due to the bottom effects mentioned, a static probe cali- 
bration method was used: With the probe in place, the water level 
in the tank was varied and the output yeas measured on the 
probe side of the blocking capacitor. (It was verified by use of 
an AN/URM-26 AM modulation test set that there was no measurable 
loss in amplitude across the blocking capacitor in the frequency 


range of interest.) 
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The Wave Height Measuring System 


Figure 10 
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The system for processing wave information is shown in 
Figure 11. The wave information, contained in the envelope of the 
wave probe signal, was first amplified; this information was then 
recorded on an FM tape recorder and simultaneously sent to a signal 
correlator. The output of the correlator was graphically recorded 
emecthie XY plotter. 


ine Sound Source wsed was a BaTiO. cylinder, 4.5 em sin 


; 
miameter and 10.5 cm tall. The receiver was a 1/8 inch BaTi0., 
probe. The sound processing system is shown in Figure 12. The 


received signal was pre-amplified in a low-noise 40 db amplifier, 
Pemapass filtered at carrier + 2.5 kHz and then amplifier further. 
The amplitude modulation of the acoustic signal was analyzed in 
like manner with the same detector network used in the wave probe 
system. This information was then recorded in the same manner as 
the wave height data. : 

The wave and acoustic amplitude fluctuation data were analyzed 
for frequency content. The system used is shown in Figure 13. The 
Original tapes, recorded in the systems previously described, were 
replayed 100 times faster and recorded on tape loops with a 


standard audio recorder. The tape loops were then replayed through 


eeo.l octave filter as shown. 


a PROCEDURE 
The experiment was conducted in water of 8.8 cm depth. A 
strong normal mode was found at 14.1 kHz, a frequency well below 


Eutcorf of the n3 = 2 family of modes. The sound field was probed 
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in order to identify the modes nodel pressures other than zero 
were found in the tank, indicating excitation of more than one 
mode or more likely imperfections in the boundaries. The cri- 
terion for a strong mode was taken to be the average standing wave 
ratio, which in this case was approximately 60. 

The mode was identified by measuring the lengths of the 
standing waves and estimating the wave numbers. The mode with 


n. = 28, n 


1 = 9, n, = 1 appeared to be excited. This mode has a 


3 
Pebculated frequency of 13.6 kHz (4.25 percent less than the 


e 


observed frequency). The vertical sound protLiule Verified, that 
the ng = 1 family of modes was being excited (Figure 14) as 
assumed in Section II. 

Once the desired mode was excited and identified, three re- 
ceiver stations were selected at strong antinodes of the hori- 
zontal wave pattern (Figure 15). The procedure now described was 
followed at each receiver position. 

The acoustic probe was set at a depth of 4.4 cm in still water, 
and the sound carrier amplitude measured. The surface was then 
disturbed with a preselected wind condition. The following ob- 
servations were then made of the envelope of the modulated sound 
carrier: 

1. A 500 second record was made of the mean square 
amplitude of the modulation from the output of the 


correlator on the XY plotter. 
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2. A 20 minute FM tape recording was made of the 
acoustic modulation. 

3. The autocorrelation function was computed at two 
different times by the correlator and recorded 
bi the Soap Letter. 

The wind generator was then stopped and the surface allowed 
to calm. The carrier amplitude was again measured to verify that 
the amplitude of the standing wave had not changed. The acoustic 
data for six different sea bares at each receiver position were 
Soerained in this fashion. | , 

The area directly under the fans of the wind generator was 
intensely agitated relative to the rest of the sea surface. In 
order to isolate this contribution to the modulation of the signal, 
the tank was modified as shown in Figure 16. The reverse flap 
lifted the wind off the surface, and the styrofoam prevented the 
propagation of surface waves. Since the styrofoam was pressure~ 
release, its presence had a negligible effect on the standing wave 
pattern. 

With the probe in position at a designated receiver station 
and the wind machine operating at a specified setting, the «mpli- 
tude of the acoustic modulation was measured through a 0.1 octave 
filter and average modulation amplitude as a function of frequency 
was recorded. These data were recorded at positions "B" and "Cc" 


Sly; position "A was obstructed by the styrofoam buffer. 
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Once all sound data had been collected, the statistics of the 
sea surface were measured. A calibration of the wave probe was 
always conducted before a data run. With the probe in place, the 
tank was filled with water to a depth of 10 cm. The water was 
then slowly siphoned out and the DC output of the wave probe was 
recorded as a function of depth. A calibration factor was de- 
termined from this curve. Probe performance did not vary with 
position in the tank except very near the walls. 

—peotatic calibration oF the wave probe is considered valid; 
however , a dynamic check was performed by simultaneously measuring 
the height of a passing wave visually and comparing it with the . 
output of the probe. The passing waves were measured on a metric 
scale, printed on heavy photographic paper and taped to the side 
of the probe holder. (The thin paper ruler caused negligible dis- 
turbance in the passing surface waves.) The output of the probe 
was observed on an oscilloscope. 

Once the probe was calibrated, it was placed at each 
receiver position and measurements were made of the amplitudes of 
the waves generated on the surface. 

The statistics of the surface at the position of the sound 
source were also required. These data could not be taken with 
the wave probe, so the behavior of the various frequency components 
of the sea were measured as they propagated down the tank anda 
projection of the behavior of each water surface at the position 
of the acoustic source was made. This estimate was checked 


visually and appeared to be valid within experimental error. 
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Analog frequency analysis of the FM tape recording was made 
by a tape loop method. The 20 minute tape records were played 
back 100 times faster and recorded on 7 second tape loops with and 
audio-magnetic tape recorder. 

These loops were then played throuch a 0.1 octave filter with 
slow integration time. Average amplitude as a function of fre- 
quency was recorded. The tape loop output was also played through 


a true RMS voltmeter to measure Overall RMS amplitude for com- 


parison with the mean square amplitude prediction of the correlator. 


D. ANALYSIS 

The RMS wave height was computed by having the correlator set 
Hoecompute the autocorrelation function. In this mode the cor- 
relator output is 

T 
c(T) = = ic Vin (t)at (3.6) 

which equals the mean square value of the input when T approaches 
ietiaity. The correlator integration time is only 20 seconds, so 
a good estimate of the mean square value was obtained by taking 
the average of a 500 second graphical record of the correlator 
output (Figure 17). The various equipment amplification factors 
were divided out and the square root taken. 

The RMS value of the envelope amplitude of the wave probe 
Signal was computed in the Same manner. 

The average surface wave and acoustic amplitudes as functions 


of frequency were obtained as previously discussed. As the data 
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were recorded, a 100 Hz tone was recorded on another channel of 
the tape. The system amplification was checked by monitoring this 
test signal at various stages in the frequency analysis. This 
check also gave assurance that the data were not being distorted 
during the processing. 

In order that a and € be computed, the gains of the equipment 
through which the data was sent had to be divided out of the 
observations. 


Refer to Figures 11 and 12 and let 


GAIN jwp = the amplifier gain 
GAIN, app 1 the FM tape recorder 
GAIN, app 5 = the audio tape recorder. 


The average amplitude measured through the 0.1 octave filter 
is converted to the average value of € for each frequency interval 


as follows: 


> 


; > 
WAVE PROBE OUTPUT = aa ee - GAIN 
AMP iVee al TAPE 2 


The wave probe output can be expressed in centimeters of 
wave height by dividing by the DC calibration factor determined 
before each run: 


PROBE OUTPUT 
CALIBRATION FACTOR 





WAVE HEIGHT = 


The normalized fluctuation € then becomes 





_ WAVE HEIGHT _ G12 SOWA BEI TL ey 1 \/1 
; : DEPTH 
DEPTH Coe ee eee Gale FACTOR] |DE i 
(357) 
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FILTER OUTPUT 1 | 


f 
Ben eee SACTOR « DEPTH GAIN. GATN.. 4° 
GAIN ap CAL FACTOR DEP EH (GAT, 1 GAIN. jpp . 


The sound amplitude fluctuation a at each frequency of 
interest is 


. ENVELOPE AMPLITUDE AT THAT FREQUENCY 
CARR TER <AMPLD TEUDE 


Since the sound information was processed with the same system 


as the surface data, a was computed as follows: 


_ AVERAGE AMPLITUDE : 1 (3.8) 
T e e 
RORTN yp) (CARR nao) (CAINS apy: (ech a . 


Notice that no calibration factor is required in the calculation 


Or Cs 


il 


rey .Ga eee J 
= TAe AAP Eee 


In the ratio a/e€, the factor vanishes; 


therefore, all data were treated in arbitrary units with this 
Hretor included. 

The carrier amplitude by which the envelope amplitude was 
divided was an average of *the amplitude measured before and after 
each run. If, however, the difference in these amplitudes was 
greater than 2 percent, the run was re-executed. 

The presence of high frequency sound eee eee unpre- 
dicted by the theory, and the persistence of these fluctuations 
throughout the tank strongly suggested the fans as a source. The 
acoustic amplitude fluctuation que to the agitation under the fans 
was obtained for positions "B" and "C'" and described above. For 


any given sea state, these data were very similar at the two 
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positions. Therefore, an average (fan noise) level for each sea 
state was assumed (Figure 18) for all the receiver stations. 

The other factor to be accounted for was source modulation. 
The sea state at the source was predicted as previously described.. 
The effect of the waves on the acoustic output of the source was 
determined by raising the source a given distance out of the water 
and measuring the reduction in sound amplitude at position "Bp", 
With a 10 percent reduction in the radiating surface, the sound 
pel itude at position ''B" ae Medticed by GO percent this relation 
was found to be linear enn the range of interest. Therefore, 


the amplitude modulation of the source was assumed to be 


an =" 0.26 Cra (329) 


(This result was checked by shading the upper part of the 
transducer the same amount with 3/8 inch neoprene foam rubber and 
comparing the results, which agreed within deer iene ae 
The amplitude fluctuation due to the local surface effect was then 


2 


computed with Eq. 3.5: 


2 2 2 
a =e, eucee. C- : (35) 


As a check on the validity of the data, several quantities 
were measured by two methods and the results compared. The tape 
loops, for instance, were played through the true RMS voltmeter in 
order to measure the RMS amplitude fluctuation and wave height. 
The results, adjusted for system gain, agreed to within 10 percent 


with the RMS values computed through the correlator. 
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The frequency distribution of the modulation data was also 
determined from the autocorrelation function plotted during each 
data run. The cosine Fourier transform of the autocorrelation 
function c(T) equals the power spectral density. The transform 
of c(T) was computed by numerical integration of 100 points taken 
from the plots described above. The resulting spectra confirmed 


the distributions obtained from the tape, loops. 


Eee ERROR ANALYSIS i: 

The major cause of error: in this experiment was the non- 
stationarity of the processes. The variation in mean Square 
amplitude (see Figure 17) gives a qualitative measure of the 
degree of non-stationarity involved. 

Measuring equipment with aera long integration times 
to overcome this problem was not available. Therefore, long 
"effective'’ integration times were achieved by visual integration 
of the meter readings or .graphical averaging as previously dis- 
cussed. Study of Figure 17 suggests that an integration time of 
500 seconds is sufficient for our purposes. 

In the eveaioncy analysis, 15 minutes of FM tape record was 
compressed onto a 7 second tape loop (see page 33). Visually 
averaging the meter reading on the 0.1 octave filter, provided an 
integration time in excess of 500 seconds. The estimated error in 
the value assigned to acoustic amplitude fluctuation (a,) ox sur - 
mace height fluctuation (€) is + 10 percent. (This estimate was 
determined by comparing the values estimated by two independent 


observers. ) 
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The tape loop method was not used in obtaining the acoustic 
amplitude fluctuation due to the fans; in this case the filter was 
placed directly on the output of the filtering and detecting 
network (page 40). Averaging two sets of observations together 
@t each wind setting, minimized the error in the fan contribution 
(a-) to an estimated 15 percent, so that the more elaborate tape 
loop analysis was not deemed necessary. 

Recall that the acoustic amplitude fluctuation term a 
resulted from the source modulation, which was extracted from 
the € curves. The Boece exror in © is + 10 percent. “Adding 
a + 5 percent error for the graphical extrapolation yielding the 
surface at the source position, results in a total error in the 
source modulation tern, aes of + 15 percent. 

The RMS amplitude values computed from the graphical output 
of the correlator are considered reliable. In order to test the 
validity of the average amplitudes obtained from the tape loops, 
RMS amplitudes computed in this manner were compared to those 
from the correlator. The values were consistently within 
10 percent.. 


mine error in Q,: a, @ and € were assumed statistically in- 


re 
dependent at each receiver position. Therefore the error in the 


value of a/€ at each frequency was computed using standard RMS 


error analysis. 
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IV. RESULTS AND CONCLUSIONS 


The exper imental wesults consist Of eidhtecnescets Of dare 
measurements were taken at three different receiver positions for 
six different winds. The six winds generated can be grouped in 
mers: (1,2), (3,4) and (5,6). Hach pair of winds was generated 
with the same settings of manifold height, angle of attack, and 
Seifice size of the wind machine. The even-numbered winds were 
generated with five fans; the odd-numbered ones, with three. The 
receiver positions were previously shown in Figure 15. 

The frequency distribution of the measured acoustic-amplitude 


fluctuation (a); the fluctuation induced by the fans (a and the 


y 
amplitude fluctuation due to source modulation (a.) are plotted 
for each wind setting at a given receiver position (see Appendix B). 
These graphs reveal certain general features: 
1. The frequency Sheen on of the measured acoustic 
GIVE tua talon (a. ) closely resembles that of the local 
surface fluctuation (€). In all cases each curve computed 
for a closely resembles the curve for the napa ileme a. 
2. At each receiver position, the frequency spectra for 
both) “Band "Care Similar for each wind pair. The 
stronger winds produce greater amplitudes as would 
be expected. 
Ss. the curves cf aq, are Similar at all receiver positions 
for any given wind setting. The difference in values 


of a, LOn pOSTELOns. 5. Andee Gals sViemy, sina lie 


pail 








4, The higher frequency components of € appear to die out 
as the surface waves progress down the tank. (See 
Figure 22.) The corresponding values of a. takiont 
much more slowly. 

5. In computing the values of the acoustic amplitude 


fluctuation resulting from the local surface 
a 2 2 2 
a = \(0.) a (eye (ol. (Bm) 


we see that a, becomes significant at the higher 
frequencies while a. 1s Significant at the lower 
f1;equencies. 
6. There appears to be evidence of generation in a, of 
second harmonics of frequencies in the region corre~ 
Sponding to resonance. 

The computed values of a/€ (using Eq. 2.13) are plotted as 
functions of surface-wave frequency for each wind condition at 
each receiver position. “(Appendix C.) The results are summarized 
in Figure 19 where the curves are collected on single graphs 
Mecordging to wind condition. 

In the region below the predicted ''y esonance" the measured 
values of a/€ compared favorably with the theoretical prediction. 
Very little evidence for a resonance effect exists in the results 
of any run. Above the predicted resonance the observed values of 
Myo are consistently greater than the theoretical values. The 
absence of any peaking of a/€ in the vicinity of resonance cannot 


be explained. 


ae 








Several important factors could be contributing to the high- 
frequency behavior mentioned. Recall that. the term ae 1S eX- 
tremely important in determining the acoustic amplitude fluctu- 
ation ain this region. [t is suspected that perhaps the entire 
observed acoustic amplitude fluctuation a. above 6 Hz can be 
attributed to the surface agitation under the fans. 

In the design of the fan isolation experiment, it was assumed 
that the surface agitation directly under the fans was responsible 
for the additional Preece. ume rarer at the receiver. The 
styrofoam buffer was therefore placed under the orifice flap. 
This, however, suppressed large surface waves generated in this . 
region. It is conceivable that these waves contribute to the 
"distant sources", but this possibility was not investigated at 
the time. Recall also that analyses of the fan-generated noise 


were not conducted with the more accurate tape-loop method 


because of time and equipment limitations. Therefore, the de 


»s 


term may be larger than recorded. This hypothesis would have 
been investigated further if time had been available. 

Since accuracy in a becomes critical when computing the 
mutao g/€ for the higher frequencies, increasingly Ronetors 
behavior in a. in this frequency range at distances further 
removed from the fans cannot be considered as unequivocably con- 
meadicting the theory. 

The discrepancy between experiment and theory for values of 
fe at the higher frequencies could also be attributed to 


violations of the assumptions in the theory. For instance, 
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Figure 20 shows that the surface fluctuations were functions of 
both x and y, rather than just of x as was’assumed. Investigation 
of a more realistic expression for F(x,y,t) would require extensive 
auto- and cross-correlation on measurements of the generated sur- 
mace fluctuation. 

It is also assumed in the theory that the properties of the 
surface fluctuations are uniform in the vicinity of the receiver. 
Figure 21 shows the distribution of RMS wave heights for a typical 
wind over the area of the Surface in the vicinity of the receiver 
stations. This figure shows that the waves are dying out as they 
propagate down the tank. Figure 22 shows that the frequency 
spectrum for the surface waves above the frequency corresponding 
to the "resonance" is also changing. The effect of the inhomo- 
geneity of the surface could not be determined but is probably 
mOSt severe at position "A", 

A basic assumption in the theory is that the mechanism for the 
Meeustic fluctuation is the local surface state. In the investi- 
gation of the validity of this assumption, the following points 
merit study: 

1. With the normal mode excited and the water surface 
calm, the surface in the vicinity of the source was 
disturbed by quadradigital Saee tea: An immediate 


amplitude fluctuation was observed at the distant 





Trickling four fingers in the water, making small capillary 
waves. 
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The Tank Surface Disturbed 
by a Typical Wind 


Figure 20 
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receiver (station "B'). The surface above the receiver 
was disturbed in the same manner and the acoustic signal 
was again immediately perturbed. However, when the water 
surface was disturbed at a point removed from both source 
and receiver, an effect on the received acoustic signal 
was observed only after the surface disturbance had 
propagated into the vicinity of the souree, receiver. 
and region in between. 

2. the treatment of the extreme agitation under the fans 
as a distant source (page 14) contributing to observed 
acoustic-amplitude fluctuation appears valid from the - 
results of the fan isolation experiment, but this would 
appear to be inconsistent with the result just discussed. 
Recall, however, that the violent agitation of the sur- 
face under the fans is a clear violation of the assumption 
that ©€F(x,y,t) “be small. 

The first point above appears to validate the assumption that 
the driving mechanism for the acoustic fluctuation at the receiver 
is the local surface. The fan isolation experiment: (point two) 
showed that severe surface disturbances at a remote eee could 
cause amplitude fluctuations under a smooth surface. It is plau- 
Sible that the distinction between these two cases results from 
the difference in magnitude of €F(x,y,t) in the two regions. 

The stringent assumptions of the simple theory were all clearly 


violated to some degree during the conduct of the experiment. It 
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is concluded, however, that the measured values of a/€ below the 
frequency of predicted "resonance" are sufficiently close to the 
theoretical values to claim consistency. between experiment and 


theory in this range. 
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APPENDIX A 
Computer Program for a/€ Calculation 


TIPLICIDZ REAL *Sieams KP.) 
REAL*8(DCOS, DCOSH, DABS, DSQRT) 
DO CO mk—l45 00 

GAM = (I-1)*.01 
K1=2.*3.1415926535/16.4 

KS = O26 59557 0.6 
PI=3.1415926535 
AR1=2.*GAM*K1 

AR2=K3**2 

AR3=GAM* *2 

AR4=PI/2. 

AR 5=(1.-((AR3+AR1) /AR2) ) 
AR6=(1.-((AR3-AR1) /AR2 ) ) 

Lh Ae oer sO. Go TO 50 

AR7=DCOS( AR4*DSQRT( AR5) ) 

GO TO 60 

AR7=DCOSH( AR4*DSQRT ( DABS( ARS) ) ) 
IF(AR6O.LT.O.) GO TO 70 
AR8=DCOS( AR4*DSQRT (ARG) ) 
AR8=DCOSH( AR4*DSQRT( DABS(AR6) ) ) 
AR9=DABS(1./AR7-1./AR8) 
ABE=PI*ARQ/4. 

PRINT, GAM, ABE 

CONTINUE 

STOP 

END 
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APPENDIX B 
Observed Data 
The eighteen sets of data described on page 51 are plotted as 
functions of surface wave frequency (§/2m). 
= the measured acoustic amplitude fluctuation 
€ = the normalized surface wave fluctuation in the 
VICINIty of the receiver 
a. = the acoustic amplitude fluctuation generated by the 
severe agitation under the fans 
a = the acoustic amplitude fluctuation due to source 
modulation 


Representative error flags for these data are presented in 


Figure Bl. 
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APPENDIX C 


Measured Values of @/€ 


The measured values of a/€ at each receiver position for 
each wind condition are plotted as functions of surface wave 


frequency (f/27). Standard RMS error flags are presented. 
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